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Abstract. We study in detail, solving the Bogoliubov-de Gennes equations, the 
magnetic field, supercurrent and order parameter profiles originated by a solenoid 
or magnetic whisker inserted in a type-II superconductor. We consider solutions of 
different vorticities, n, in the various cases. The results confirm the connection between 
the vorticity, the internal currents and the boundstates in a self-consistent way. The 
number of boundstates is given by the vorticity of the phase of the gap function as in 
the case with no external solenoid. In the limiting case of an infinitely thin solenoid, 
like a Dirac string, the solution is qualitatively different. The quasiparticlc spectrum 
and wave functions are a function of n — n ext , where n ex t is the vorticity of the solenoid. 
The flux is in all cases determined by the vorticity of the gap function. 



1. Introduction 

The effect of a magnetic field on a superconductor has attracted interest for a long time. 
At small fields the superconductor is rigid to the external field but for a sufficiently strong 
field, or a high enough electric current flowing through the material, superconductivity 
is destroyed. In type-II superconductors there is a low critical field, H cl , above which 
the field lines penetrate the superconductor in the form of quantized vortex lines. As the 
field increases, the vortex line density also increases until the vortex cores overlap and the 
system becomes non-superconducting at the high critical field, H C 2. The external field is 
shielded by the appearance of compensating currents that perfectly cancel the external 
field beyond a penetration length, A. Another length which is important is the coherence 
length, (, which measures the range of the establishment of the superconducting order 
parameter into the superconducting region. Since many important superconductors, 
like the high-T c materials, are strong type-II superconductors, the presence of external 
magnetic fields implies proliferation of vortices and, therefore, their influence has been 
thoroughly studied. 

In this work we consider the effects of a vortex induced by the insertion of a 
long solenoid or magnetic whisker of very small width in the superconductor. Due 
to recent technological advances it is possible to construct magnetic systems of nano 
size of the order of the coherence length or penetration length. We consider that the 
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inserted foreign system is long enough so that the field lines penetrate the material 
(eventually through the effect of the vector potential) but that the field lines return far 
from the superconducting film, such that no antivortices are created inside the material. 
Solving the Bogoliubov-de Gennes equations we study in detail the energy spectra and 
its consequences on the physical properties due to the presence of the external solenoid, 
paying special attention to the magnetic shielding due to the Meissner effect. 

The problem of the quasi-particle states due to the presence of a vortex in a s-wave 
superconductor was solved long ago both analytically pQ and numerically PJ. There 
are bound states localized in the vicinity of the vortex location and a continuum of 
delocalized states. The case of a rf-wave superconductor lead to some controversy but it 
was established that the states are delocalized, consistently with a gapless spectrum [3]. 
Classifying the states in terms of the angular momentum around the vortex line, allowed 
to determine that there is a branch of boundstates, one for each angular momentum value 

The results were obtained looking for an order parameter of the form A = A e~ inip 
where <p is the polar angle and n fixes the vorticity, chosen originally as n — 1. The core 
states are coherent superpositions of particle and hole states and interpreted as being 
the result of constructive interference of multiple Andreev scattering from the spatial 
variation of the order parameter Also, it was shown that the main contribution to 
the supercurrent is originated in these states. At very low temperatures the quantum 
nature of the bound states gains importance and oscillations in the various quantities 
are quite pronounced [5]- The quantum limit is obtained when the thermal width is 
smaller than the level spacing and is reached when T/T c < l/(kp(o), where T c is the 
critical temperature, kp is the Fermi momentum and Co — Vp/A is the coherence 
length. In the typical type-II layered superconductor NbSez the critical temperature is 
T c = 7.2K and kp(o ~ 70. The quantum limit is reached for rather low temperatures 
of the order of 50mK, but for the high temperature superconductors, where the critical 
temperature is high and the coherence length is small, the quantum limit is reached for 
temperatures of the order of 10K. In this limit, Friedel-like oscillations are found in 
contrast to a Ginzburg-Landau description. In particular the Kramer-Pesch effect [6 j 
where the vortex core size increases with temperature, is not explainable in terms of a 
description in terms of normal electrons. The oscillations in the various quantities are 
observed at sufficiently low temperatures, irrespective of the size of kp( . 

Even though as the field increases it is energetically favorable that many single 
vortices appear, as compared to fewer vortices with higher units of flux j^j, the possibility 
of multi-flux vortices (or giant vortices) has been considered in particular systems and 
situations. In this case it has been predicted that there are n branches of boundstates 
(where n is the vorticity of the flux line) [S] and confirmed by several authors j^j. 
If formed, the doubly quantized vortex is metastable against dissociation into singly- 
quantized vortices. In this case a counter-circulating current appears in the core [I] 
which is due to a bound state that appears below the Fermi level. At high temperatures 
this state is depopulated, the current reversal disappears and leads to a structure similar 
to the Ginzburg-Landau description. The appearance of n discrete branches of states 
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for a vortex with vorticity n should be seen in Scanning Tunneling Microscopy (STM) 
measurements in the form of n rows of peaks as a function of the distance from the 
vortex core [Oj. Oscillations in the current due to the multiple vortices were observed 
at low temperatures in a form that is qualitatively different from results obtained from 
a Ginzburg-Landau description ^0]- These oscillations are also observed in the order 
parameter. 

Multiply quantized vortices have been observed in type-I superconductor films 
at high magnetic fields ^Tj. In high temperature superconductors having columnar 
or large pointlike defects acting as pinning centers, multiply quantized vortices may 
appear as well. Giant-vortex states have also been observed in superconductors with 
strong geometrical constraints, like in triangles or squares of sizes of the order of the 
fxm, using a scanning SQUID microscope ^2j. The possibility of different vorticities 
leads to interesting features like the Little Parks effect where oscillations of the critical 
temperature as the external field changes occur due to transitions between states with 
different vorticities [T3] . 

In general, we may have a mixture of vortices with different winding numbers. 
Recently, in a NbSe2 traditional superconductor, where normal metal islands of gold 
are inserted, a coexistence of strongly interacting multiquanta vortices distributed in 
a lattice with interstitial single vortices, has been observed (T^j. The case of a single 
normal dot inserted in a superconductor was also studied [TH]. The dot is quite small 
of a subnano meter scale. It was found that the energy spectrum strongly depends on 
the number of flux quanta penetrating the dot and it was shown that the number of 
branches corresponds to the number of flux quanta. To guarantee that any vortices in 
the system would be in the dot, fields smaller than H c i were used. It was shown that by 
increasing the dot size the vorticity of the dot could be increased and therefore different 
states could be studied directly. 

The vortex lines in general appear due to the application of an external magnetic 
field typically homogeneous. However, we can as well consider the presence of magnetic 
field lines that are due to a solenoid or a magnetic rod inserted in the superconductor. 
Actually, the magnetic field lines do not need to penetrate the superconductor itself, 
since what really matters is the vector potential. It is the vector potential that appears 
in the Hamiltonian of the system in the presence of a magnetic field [Zj, as is well known. 
This has been emphasized [TH] considering a superconductor in the form of a cylindrical 
shell of internal radius R and width a in the center of which is inserted either a solenoid 
or a magnetic rod of radius r smaller than R. In these systems, considering the length 
of the cylinders very large compared to the radius, the field lines will close far from the 
supercondutor and therefore no field lines penetrate the cylinder. However, the vector 
potential due to the flux contained in the transverse section is non-zero and the field 
has the same effect on the super currents. 

We may as well consider other situations in which there is an interplay between 
magnetic systems and superconductors. Recently the situation of ferromagnet- 
superconductor hybrid systems has been reviewed [T7]. On one hand it has been 
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proposed that we can manipulate spin and charge in magnetic semi-conductors using 
superconducting vortices, with applications in spintronics ^H|- On the other hand 
we may consider the effect of the magnetic system on the superconductor. The 
magnetic film is separated from the superconducting layer to avoid suppression of 
superconductivity by the proximity effect. This is attained introducing between the 
two films a thin insulating layer. 

A possibility to insert a magnetic field in the superconductor is through the field 
originated by a magnetic dot. Depending on the strength of the magnetization of the 
dot we may have single or multiple vortices. Using a London theory it has been found 
that giant vortices occur when the dot size is small enough ^H] with a size of the 
order of 4.5£. Otherwise the energetically favorable situation is the presence of single 
quantized vortices [19J. Considering an array of magnetic dots these originate vortices 
in the superconductor which are preferably bound to the magnetic dots in a way that is 
more advantageous than due to the usual defect pinning centers [20]. In zero field the 
dots, which are far apart, are not coupled and are oriented randomly. Therefore any 
other vortices that may not be bound to the dots or that appear as fluctuations will feel 
the presence of a random magnetization and will be in a resistive state. However, if a 
field is applied the dots align and if field cooled the film may be superconducting. The 
ordered commensurate state is then favored. The pinning turns out to be more effective 
than with non-magnetic dots. Also, the critical current has the unusual property that 
it is increased with field. Therefore introducing magnetic nanoparticles or nanorods 
one obtains what are called frozen flux superconductors [2H 122] • In these systems a 
thermodynamically stable state of frozen flux lines is obtained showing that the magnetic 
dots are more effective than other pinning centers. This is important to control the 
maximum currents that may flow in the superconductor. Additional flux lines created 
by an applied magnetic field need to overcome big energy barriers in order to move. As 
mentionedabove larger systems were also considered such as two films, one magnetic and 
the other superconductor, close by [23]. In this close vicinity the magnetic field created 
by the supercurrents interacts with the magnetic subsystem. The interplay between the 
two systems leads to interesting magnetic structures like the cryptoferromagnetic state 
|24"j and commensurability effects lead to increase of critical current |25j . 

Since the total flux created by the magnetic dot is zero if the superconducting film 
is large enough, one may expect the presence of vortices and anti-vortices, due to the 
dipole field of the dot. In a type-II superconductor the interaction between two vortices is 
repulsive and the interaction between a vortex and an anti- vortex is attractive. In studies 
close to the critical line it has been found that the antivortices may coexist with the 
vortices but away from the vicinity of the critical line they disappear. However, in type-I 
superconductors the interaction between a vortex and an anti- vortex is repulsive (and the 
interaction between two vortices is attractive). Therefore it is possible that structures 
with vortices and antivortices may appear in type-I systems. Also, small systems show 
confinement effects on the vortices: due to the strong increase of the kinetic energy term 
near the frontier, it is favorable for the system to be in the superconducting state close 
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to the borders. Therefore it is not favorable for a vortex to approach the borders of the 
system. As a consequence they tend to be close together near the center, particularly if 
the system size is very small. Considering a type-I system with a triangular symmetry 
it has been shown that indeed stable vortex-antivortex structures are possible [21]. A 
rich variety of vortex structures has been studied recently fI7\ . 

All of these results show the importance of the study of complex structures of 
vortices and their interplay with magnetic structures. Recently it was stressed that the 
important characteristic that determines the boundstates is the winding of the phase 
[2E] • The detailed form of the order parameter in the vicinity of the vortex core is not so 
relevant. Performing a non- self-consistent study of the spectrum it was found that the 
suppression of the gap function has a minor role. The important feature is the winding 
of the phase. The supercurrent acts in non-symmetrical way on the particle and hole 
parts of the quasiparticles. It tends to decrease the angular momentum of the particle 
part and to increase the angular momentum of the hole part. 

In this work we consider the same problem from a different point of view. We 
perform a self-consistent solution of the influence of a very long solenoid on the properties 
of the system. We consider different possible vorticities and study how the system 
responds to the external perturbation. The internal field must adjust itself to the 
solution chosen according to the external field exerted by the solenoid. The total 
magnetic flux is fully determined by the choice of the angular momentum of the gap 
function and the value of n determines the vorticity of the vortex solution. This may 
be a single or a multiple vortex. Depending on the relation of the value of n and 
the value of the external field, the internal currents will create fields that compensate, 
undercompensate or overcompensate the external field. The various situations lead to 
different spectral structures depending on the width of the solenoid. These in turn 
originate different structures for the internal field and supercurrents generated. The 
limit of a very thin solenoid is qualitatively different. It is shown that when the vorticity 
chosen equals the unit of external flux the currents generated vanish and no bound states 
appear. Otherwise the currents may be positive, and branches of boundstates with 
positive energies appear or the currents are negative and branches of negative energies 
appear. These results confirm the link between the boundstates and the internal currents 
in a self-consistent way. In the case of a finite width solenoid the number of boundstates 
equals the vorticity of the gap function and is insensitive to the external field. 

2. Method 

2.1. Bogoliubov-de Gennes equations 

Consider the Bogoliubov-de Gennes equations (BdG) 
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where U(r) is an external potential, A(r) is the vector potential and where we consider 
s-wave pairing, for simplicity. A(r) is the pairing function given by 

Mr) =9 E «i(rK(r) [l-2/(^)] (2) 

Here f(Ei) is the Fermi-Dirac distribution. The vector potential is given by Maxwell's 
equations 

4:71 

VxB = VxVxA = — 3 total (3) 

c 

which, in the Coulomb gauge ( V.A = ), is given by 

4-7T 

V 2 A = -^Jtotal (4) 

The current density originated in the supercurrents is obtained self-consistently by 

Ui(r) 



eh 



V - f A(r) 

he 



+ [(1 - /(£yh(r) 



v*(r) - c.c. (5) 



2.2. Diagonalization of the Hamiltonian 

We assume no dependence along the axis of the vortex line (z-axis) and cylindrical 
symmetry both of the superconductor and of the potential U(r). Let us take the order 
parameter in the form 

A(r) = A(p)e"^ (6) 



This form describes a magnetic flux equal to n flux quanta ( $ = n$o = n 
The wave functions Wj and are expanded in a way similar to ref. [2] 



he 
2c 



^( r ) = E^^-i/2^- 1/2) ^ (7) 
<*) = Y.<i^-V^- 1/2 ^ (8) 



where the basis functions are 



= ^7 7- — r^m ( a jm — ) (9) 



RJm+l{.Q.jm) V -R/ 

The system is placed in a cylinder of radius i£. Given the azymuthal symmetry of 
the system, neither A(p) nor A depend on (p. Therefore the Hamiltonian may be 
diagonalized separately for each value of the angular momentum //. The functions J m 
are the spherical Bessel functions and <x,- m is the j th zero of the Bessel function of order 
m. The set of values of the angular momentum is given by fi = ±(2/ + l)/2 where 
I = 0, 1, 2, ■ • -. The terms (fi — 1/2) and (jj, — 1/2 + n) may be rewritten in a more 
symmetrical way like (// — n/2) and (jj + n/2) if jj! is half-odd integer if n is odd and 
fi' is integer if n is even. 
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For each eigenvalue E iy we have a single value of p and it is enough to diagonalize 
the matrix, defined in the subspace of the zeros of the Bessel function, 



where 
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+ ^n-E F d2) 

- [ R A (p) 

h = / (t>j,^-i/2(p)^ E —(t> j '^-i/2(p)pdp (13) 

/■-R A (p) 

A + = y o 0i,M-l/2+n(p)^^0i',^-l/2+n(p)prfp (14) 
J 2~ = / 0i,M-l/2(p)^(p) 2 0i',^-l/2(p)p^P (15) 

f R 

U = / 0i,M-l/2+n(p)^(p) 2 0i'^-l/2+n(p)P^P (16) 



Also we have 



0j>-i/2(p)|A(p)|^v iM _i /2+n (p)prfp (17) 
It is important to note that the symmetry of the BdG equations 

Ui(r) -> t£(r) (18) 

Vi (r) - -<(r) (19) 

^ -> (20) 

allows to reduce the solution to the positive values of p. We obtain the eigenvectors and 
eigenvalues for negative values of p using the above symmetry. 
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2. 3. Calculation of the vector potential 

Consider a general situation where the total vector potential is given by the sum of 
an external potential and the internal vector potential originated on the supercurrents. 
Then we can write that 



A = A ext + a 

where a is the internal gauge potential. Then naturally we can write that 
$ = $ a + $ ext 
B =VxA = Vxo + B ext 
V 2 A = V 2 a + V 2 A ext = V 2 a 



(21) 



(22) 



since V 2 A ext = 0, except at the region where the external currents are non-zero. 
Therefore we have to solve the equation 



V 2 a = -^J 



(23) 



Due to the quantization condition the total flux is given by 

$ = $ a + <$> ext = n$ (24) 

Therefore we will be considering situations where the supercurrents will generate internal 
magnetic fluxes such that 

$ a = (n - n ext )®o 

where & ext = n ext $o (we take real parameter). 

Let us then start from the equation 

1 d ( da u 
pdp 



(25) 



P- 



47T 

J u-. 



dp J p" c 1 
Defining a v = F(p)/p, Poisson's equation reduces to 

d 2 F ldF 4tt 
dp 2 p dp c * 
Since the current is given by 



(26) 



(27) 



J v {p,z) = - £ f(Ei)\ui(v) 



+ [l-f(E t )}\ Vl (r 



we can make the decomposition 
An 



i,/i>0 
2 



p-1/2 



p-l/2 + n 



Yc A ^ z) 



— J^p = K(p)+P(p)F(p) 



(28) 
(29) 
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with 



4-7T 

k(p) = -— E/MklV-i/2) 

mc i 

- [l-/(^)]|^| 2 (/i-l/2 + n) 

+ ^p a ?(p) E {/Wkl 2 + [i - /WIN 2 } 



and 



mc he 

Therefore we get 

d 2 F 1 dF 



dp 2 p dp 
Discretizing this equation 

d 2 F F i+1 -2F i + F i - l 



dp 2 a 



2 



dF 



i+i — -fj-i 



(30) 



P(p)--^MT,f( E i)\u i \ 2 + [l-f(E i )]\v i \ 2 (31) 



= K(p) + (3(p)F(p) (32) 



dp 2a 

where a = R/ (Nr + 1), and N r + 1 is the number of points. We get 

F ' +1 " + ^ - Fi+ \~ Fl ~ l — Ki + m (33) 
a z 2pa 

The boundary conditions are such that F(0) = 0(0^ does not diverge at the origin) 
and ^\ p=R = ( B = outside). In the first case it is enough to take Fq = 0. In the 
second case we have F N+ i — F N = 0. In the outside boundary, i = N r , we have 
-F N + F N ^ F N - F/y_i 
a 1 2a(R — a) 

The system is therefore reduced to a tridiagonal system of equations. 

The equations are solved self-consistently choosing appropriately the first iteration. 
The solution converges after a few iterations usually less than 10 iterations. 



3. Results 

The solution of the BdG equations gives full information about the superconductor 
within BCS theory. The equations are solved in atomic units where m — h — e — 1, 
c = - ~ 137. We consider two sets of parameters. One set (designated set A) 
corresponds to an extreme case in the quantum limit and another set (designated set B) 
corresponds approximately to the parameters suitable for the traditional superconductor 
NbSe 2 - In the first case we consider parameters such that g = 0.8, E F = 0.5, uo D = 0.25, 
R = 80 — 250, n r = 1001, rij = 200 and = 600. According to these parameters we 
have that Co = vf/^o ~ 7.69, kp(o ~ 7.7. The critical temperature is of the order of 
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Figure 1. Comparison of the gap function and the magnetic field for the two sets of 
parameters A and B, defined in the text, at a very low temperature. We take n = 1 
and no external solenoid (n ext = 0). The vertical axis is scaled to the bulk value of the 
gap function and the value of the field at the origin, respectively, and the horizontal 
axis is scaled to the coherence length and the penetration length, respectively, of each 
set of parameters. 

0.1 (all of these numbers are in atomic units). The second set of parameters is given 
by g = 0.31, E F = 37.2meV, u D = 3meV, R = 10000, n r = 1001, rij = 200 and 
n M = 600. According to these parameters we have that Co = vf/A ~ 876, kp( ~ 65. 
The critical temperature is of the order of 8K. (Recall that one atomic unit of distance 
is ~ 0.5 Angstrom, and that one atomic unit of energy is ~ 27eV). In spite of the large 
difference between the two sets of parameters we will see that many of the results are 
smilar. In the first case we consider temperatures that are rather small of the order 
of T/T c ~ 0.001 which corresponds to the quantum limit regime (in the case of NbSe2 
this regime sets at a temperature of the order of 50mi^). The scale of the energies is 
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such that in the first case A ~ 3.53eV^ while in the second set it is of the order of 
A ~ 1.2meV. 

3.1. A ext (r) = 

We begin by reviewing the case where there is no external potential U(r) and no explicit 
A ext (r), in addition to a vortex line characterized by the vorticity n. In the first iteration 
we choose A(r) = A tanh £ and A(r) = A ext (r) = 0. This is the usual way to study the 
presence of a vortex line induced in the superconductor by the action of an external field, 
probably uniform. The energy spectra for different vorticities show the characteristic 
energy gap considering s-wave symmetry. For each angular momentum value there is a 
set of energy eigenvalues (in number given by the basis set of Bessel function zeros) with 
positive and negative values. The negative values give information about the positive 
energy values for the negative angular momenta, as explained above. The n = solution 
corresponds to the absence of a vortex and is therefore the result expected in the bulk 
of the superconductor away from the vicinity of any vortex. For n > 1 the solution 
corresponds to a vortex with increasing vorticity. In these cases a set of bound states 
appears in addition to the continuum. The number of states for each angular momentum 
value is given by the vorticity. For n — 1 the boundstates are all positive (the Caroli, 
de Gennes and Matricon states 0) but for n = 2, and n = 3 some of the boundstates 
have negative energies. 

The range over which the gap function reaches its bulk value defines the coherence 
length. For n = 1 the gap function is always positive but for n = 2 and n — 3 it 
has a node. This is associated with the appearance of negative energy boundstates. 
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The oscillations in the various curves are eliminated if the temperature is increased and 
in general the coherence length increases with the temperature (Kramer- Pesch effect). 
Also the coherence length increases with the vorticity. 

The magnetic field is shielded in the superconductor. In the case of a n = 1 vortex 
the magnetic field decreases monotonically evolving in the bulk to an exponential decay, 
which defines the penetration length, A. As the vorticity increases the magnetic field 
profile is no longer monotonic but there is a maximum at a finite distance from the 
vortex. Also the penetration length increases with the vorticity, as does the coherence 
length. In Fig. 1 we compare the normalized gap functions and the normalized magnetic 
field for the two sets of parameters considered. In the first case we estimate a coherence 
length of the order of £ ~ 2 and in the second case £ ~ 100. From the decay of the 
magnetic field we estimate that in the first case A ~ 43 while in the second case A ~ 790. 
This implies ratios of the order of X/( ~ 23 and X/( ~ 8, respectively. 

In order for the magnetic field to be shielded an internal field is generated in the 
superconductor to compensate the external field. This field is generated by supercurrents 
that appear in the superconductor due to the motion of the Cooper pairs. At the origin 
there is no current since the total field equals the external field and no shielding takes 
place. Considering first the n = 1 case the current increases to compensate the external 
field and then decreases exponentially. The cases of higher vorticity are more complex 
[3]. The current is negative near the origin, then becomes positive and follows the same 
trend as for the singly quantized vortex. The increased winding of the phase around the 
vortex induces in the immediate vicinity of the vortex core currents that have opposite 
circulation. In Fig. 2 we show, for completeness, the magnetic field and the flux for the 
cases of n — 1 and n = 2 for NbSe2 (parameters B). 



3.2. External field line: solenoid of negligible width 

Consider now that in addition or instead of the uniform external field we insert in the 
superconductor a solenoid of negligible width. A current may be applied and the solenoid 
(assumed infinitely long) creates a magnetic field line of negligible width inside it and a 
zero field outside. The solution of the BdG equations depends both on the vorticity of 
the gap function and on the total field through the vector potential A = a + A ext . The 
external vector potential due to the thin solenoid is given in this limitting case by 

ALt = ^ (35) 

However, in the numerical solution we replace p — > p, where p = V p 2 + 5' 2 is introduced 
to regularize the vector potential at the origin (5 is infinitesimal). In most cases the term 
which includes the vector potential may be neglected in the BdG equations because it 
is either negligible or very small and the main influence of the vector potential is in the 
expression for the current. If there is no external vector potential the internal vector 
potential (and consequently the total vector potential) vanishes at the axis (like ~ p). 
However, if the external potential has a singular form this term is important and must 
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Figure 3. Magnetic flux (internal, external and total), for n = (top panel) and for 
n = 1 (bottom panel) and different n ext in NbSe2 (set B), with T = IK. 
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be taken into account. We note that we include the vector potential term in the BdG 
equations even though its contribution is small [2|. 

When the external solenoid is inserted and a field generated, one expects the right 
solution to be the one corresponding to a vorticity matching the total field. One may 
however look for other solutions of the BdG equations which eventually will have higher 
free energies as compared to the right solution. We will in the following consider different 
possible solutions of the BdG equations even though only one is the lowest free energy. 

Consider again the BdG equations. Assuming that A(r) = A v (p)e v with 

Tic 

A<p = n ^^T p + a v ( 36 ) 

and making the general assumption of azimutal symmetry, we can write the eigenvectors 
as 

ui = g^{pV itx '- n/2)tp (37) 

Vi = h^{p)e^' +n ' 2 ^ (38) 

Then we have 

/ h 2 A d dg^j (a'-n/2) 2 2(// - n/2)e ou 

H ° Ui = ( " ^\-pTp {p -dT ) ? — 9 " s he — 7 Vj 

(fjf - n/2)n ext 

^ 9tI ' J ~ J? 9 "'* ~ ~hc~J 9fl ' J ~ W? %9>1 ' 3 < 

- Ejrg^eW-"™* (39) 
Now, grouping and simplifying the terms the expression simplifies to 

/ h 2 A d dg^j (jjf -n a /2f e a v 

= ( " 2^ ? 9 *> - ^ ~ na) ¥cJ 9 ^ 

-^<Vi)-^Vi)^'- n/2) ^ (40) 

where n a = n — n ex t. For H*Vi, the expression can be obtained from this one, by 
replacing Ui by Vi and A by —A. The first replacement is equivalent to making 

9 pi j -»• h^j,n -> -n (41) 
and the second is equivalent to 

n ex t -> -n ex t, a<p -> -a v (42) 
So, H*v is given by 



e 2 



/ h 2 A d dh^j {p! + n a /2) 2 , e a v 

H ^ = ( " 2^ 7 K ' j + {2fl + na) hcJ h ^ 

a% - E F hrt)eW+«n* (43) 



So the explicit dependence on n and n ext , comes only through n a . 
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The same choices for the vector potential and for the eigenfunctions leads to the 
current in the form 

7 e ^Vf/p \ 2 (H 1 -n a /2 e x 

+ (i - Ke^K'j C' ~ p nJ2 + Yc a «) (44) 

Also, the Maxwell equation becomes 

1 d AneTi ^ 2 /// - n a /2 e n 

+ (l-/(^))^( "^7 a/2 + ^) («) 

So neither the BdG equations, nor the Maxwell equation depends on n or n ext , separetely, 
but rather on the difference n a = n — n ext . The solution of the spectrum, eigenfunctions 
and all quantities derived from these are only a function of n a . The singular term of 
the external field line has the effect of renormalizing the states of the Bessel function 
basis if the external flux is a multiple of the quantum of flux. In particular, when 
n — n ext = n a = the system has an effective zero vorticity. However, the total flux is 
determined by n, as we will see. 

If we consider the combined transformations n a — > — n a , — > — and // — > —[/,' 
the equations are invariant. This is equivalent, as we saw before, to make the combined 
transformations n a — > — n a , a v — > — a v , and Ei — > —Ei, g — > /i and ft, — > — ry. Under 
these transformations the gap function remains invariant and the current changes sign, 
as expected. 

In Fig. 3 we plot the flux as a function of distance for two different vorticities n = 
and n — 1 for different applied fields. The total field is determined by the choice of 
vorticity of the solution: in the case of n = the total flux in the bulk vanishes and 
in the case of n = 1 the total flux in the bulk equals a quantum of flux. The external 
field is determined by the current that goes through the infinitely thin solenoid. This 
is parametrized by the flux Q ext = 0.5, 1, 1.5, 2$ . The fully self-consistent solution of 
the BdG equations yields the internal vector potential a generated by the supercurrents 
carried by the quasiparticles. The flux originated by the internal field compensates the 
external flux in order to give the correct total flux fixed by the chosen vorticity, n. As 
we can see in Fig. 3, for all values of the external field, the internal field is negative in 
order to compensate and give a zero total flux in the bulk. In the case of n — 1 the 
value of $ ea; t = 0.5$o is smaller than the total flux and the internal vector potential is 
positive. For $ ext = $ the internal field is zero and for the higher values of Q ext the 
internal field is again negative, as explained above. 

In Fig. 4 we show the energy spectra corresponding to the cases discussed in Fig. 
3 and compare with the case with no external line {n ext = 0). When the external field 
is higher than the total field and, therefore, the internal field is negative, there appear 
boundstates of negative energy in number n ext — n. When the external field is smaller, 
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Figure 4. Energy spectra for (set A) (top) different vorticities n — 0, n = l,n = 
2,n — 3, with n ext — 0; (middle) n = for various external fields 0.5,1.1.5,2. 
Note the appearance of negative energy boundstates as the field increases; (bottom) 
n = 1 for various external fields 0.5, 1.1.5, 2. Note the appearance of a positive energy 
boundstate at a field smaller than the vorticity and negative energy boundstates as 
the field increases. Also note the absence of boundstates for a field 1. 




Figure 5. Gap function for NbSe2 (set B) for (top) n = 
and for (bottom) n ex t = 1 and different vorticities n. 



1 and different values of n ex t 
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Figure 6. a) Magnetic field profiles for the various external field values for n — 0, n = 1 
and for set A. Note the opposing magnetic fields when n ex t > n changes to n ext < n. b) 
Induced supercurrents for the various external fields for the solutions with n = 1, n = 2 
for set A. Note the changes of signal close to the vortex location. 
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positive energy boundstates appear in number n — n ext . When the two numbers match 
there are no bound states, corresponding to a zero internal field. The boundstates are 
therefore associated to a non-vanishing internal field. Positive energies correspond to 
a positive internal field and negative energies to a negative internal field. When n ext 
is an integer the external field can be absorbed in the basis functions and there is no 
remaining field affecting the BdG equations (except for the internal field, which gives a 
small contribution). This is determined by n a . If n a = there is no field left (external 
or internal) and this is equivalent to a system with no vortex. In this case there are no 
boundstates, the gap function is uniform and there are no supercurrents. 

In Fig. 5 we show the influence of the external field line on the gap function. The 
results are presented for the case of NbSe 2 . When n ext = the behavior was previously 
studied. Since the system is far from the quantum limit the gap function has a smooth 
behavior. When the solenoid external field is superimposed, the shape is altered. As 
mentioned above the shape of the gap function depends essentially on the difference 
n — next- The coherence length shortens which may be possibly interpreted as a lower 
temperature. Also, the oscillations increase. 

A consistent explanation can be obtained looking at the magnetic field and 
supercurrent profiles. In Fig. 6a we compare the magnetic field profiles and in Fig. 
6b we compare the supercurrent profiles. In this last case we also consider the case of 
n — 2. As stated above, when the external field is larger than the total field, fixed by 
the vorticity of the solution, the magnetic field is negative near the vortex and tends 
to zero exponentially from negative values. When the external field n ex t = n the field 
is zero and when n ext < n the field is positive as in the case of the vortex line. This is 
clearly shown in Fig. 6. The behaviour of the induced supercurrents is more complex. 
It is also a function of n — n ext . When this difference is zero the current is zero. When 
n ~ n ext > we have a behaviour similar to the one for the vortex line with no external 
field: when n — n ext = 1 the current is positive going through a maximum and then 
decreasing, and for larger values of n — n ext the current has a node. When n — n ext < 
the current is always negative. 



3.3. External solenoid of finite width 

Consider now a solenoid of finite width. As mentioned above it could also be a magnetic 
whisker inserted in the material. The width of the solenoid is given by R s . We consider 
two cases R s — 1 and R s = 50. The first case is very similar to the thin solenoid since 
the coherence length is of the same order. However, the thicker solenoid has a width 
that is considerably larger than the coherence length and is of the same order as the 
penetration length. The external vector potential is now given by 

Kxt = e(R s - p)^ P + 9(p - R s )^0 (46) 

Also, we consider two cases. The first case is achieved considering an external 
potential U(r) that is strongly repulsive in the region occupied by the solenoid. This 
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Figure 7. Energy spectra for two solenoids with and without repulsion for different 
vorticities and external fields for set A. 
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Figure 8. Gap function for two solenoids with and without repulsion for different 
vorticities and external fields for set A. 



effectively restricts the presence of the superconducting pairs in the solenoid. A value 
of U(r) = UQ(R S — p), with Z7" = 50 is quite efficient and is such that the material 
inserted in the superconductor is opaque and the electrons can not penetrate. In the 
other case we take U = 0, and therefore the electrons are allowed to penetrate the finite 
width region where a magnetic field is inserted in the superconductor. This allows us 
to isolate the effect of the finite width region with magnetic field on the superconductor 
from the inherent space constraint due to the physical presence of the solenoid. 

Consider first the energy spectra. The case of R s = 1 is qualitatively similar to 
the vortex line and there is a set of boundstates as before characterized by n a . Even 
though the vector potential is not singular if the width is very small there is a large 
component of the vector potential and the numerical solution does not distinguish this 
case from a trully singular potential. When the thickness is larger than the coherence 
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Figure 9. Magnetic field profiles for two solenoids with and without repulsion for 
different vorticities and external fields for set A. 
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Figure 10. Flux density as a function of distance for n = 1 and external fields 
A x — 0, A x = 1 for different widths of the opaque solenoid for set A. 

length the structure of the boundstates is changed. In the case of the opaque solenoid 
the boundstates disappear altogether since the pair density near the origin vanishes, 
due to the repulsive potential. If the solenoid is transparent, the electrons can probe 
the solenoid core. However, since the magnetic field is spread through a finite region, 
the spectrum depends mainly on the vorticity and only weakly on the external flux. 
The effect is not strong enough to change the nature of the boundstates. In this case 
ther is no singular term and the spectrum is not qualitatively changed. The number of 
boundstates is therefore characterized by the vorticity of the gap function, n. This is 
shown in Fig. 8 where we compare the energy spectra for the two cases of U = 50 and 
U = 0. The effect of the solenoid on the gap function is shown in Fig. 9 where the two 
solenoids are also compared. 

The effects of the solenoid are best seen in the magnetic field profiles. In Fig. 10 we 
compare the magnetic field for the two solenoids for different cases. When the vorticity 
n < n ext the magnetic field for p > Rs is negative, to compensate the field introduced by 
the external solenoid. This does not happen when n > n ext since the external solenoid 
field does not exceed the field that would correspond to the vorticity selected. Also, 
note the difference at small distances between the opaque and the transparent solenoid. 
The field at small distances in the opaque solenoid is determined by n ext which creates a 
uniform field and determined by the supercurrents distributed for p > R s , which create 
an opposing uniform field inside the solenoid. In the case of the transparent solenoid 
the Cooper pairs can penetrate the small distance region and feel the singular nature of 
the vortex line implied by the selected form of the gap function. 

In Fig. 11 we show the flux density for the opaque solenoid for n = 1 with 
and without current flowing through the solenoid, for the two different widths. When 
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n ext = 1 the flux saturates in the vicinity of the core for the smaller solenoid. However, 
in the case of the larger solenoid the flux just follows the usual classical trend due to 
the constant value of the magnetic field inside the solenoid, assumed of infinite length. 
In the vortex line the field decreases from the core while in the solenoid it is constant. 

4. Conclusions 

In this work we have revisited the problem of a vortex in a superconductor. There 
are several reasons that took us to carry out this study. First of all there has been an 
increased interest in the interplay of magnetic and superconducting materials both from 
the point of view of the influence of the superconductivity on the magnetic materials but 
also the influence of the magnetic materials on the superconductivity. In particular, the 
vicinity of magnetic dots near the superconductor may serve as pinning centres and the 
motion of vortices in the superconductor may induce the motion of magnetic domain 
walls in magnetic materials. These properties may be of importance for the control of 
dissipation in the superconductor and the control of magnetic registers in the magnetic 
materials, respectively. The studies carried out before have shown the possibility of 
observation of giant vortices and therefore we have considered in this work different 
vorticities and studied the current and gap function profiles. 

An interesting problem to be considered is the penetration of a magnetic rod or 
solenoid in the bulk of the superconductor. We have considered this situation and 
studied the response of the superconductor as a function of the vorticity around the 
external field line. We have confirmed by a full self-consistent calculation that the 
winding of the phase of the gap function determines the properties of the system with 
the notable exception of an infinitely thin solenoid. In this case there is an effective 
vorticity, n a , which determines the nature of the spectrum. We established a detailed 
connection between the vorticity, the induced internal currents and the energy spectrum 
structure. 

An extension of this problem is to consider a superconductor in a finite slab and to 
insert two magnetic rods from opposite sides with either the same or opposite polarities. 
At the tip of each rod (or solenoid) the magnetic field lines will either leave or enter the 
supercondutor. Due to the Meissner effect these field lines will be confined to flux tubes 
in a way similar to the confinement of the chromodynamic field. Far from the solenoid 
ends the system looks very much like the problem studied here: a single infinite solenoid 
inserted in a superconductor. The behaviour near the solenoid tips is quite interesting 
for it provides an explicit realization of the confinement problem. This problem will be 
considered elsewhere. 
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